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Our interest goes to the different virial contributions to the equation of state of charged colloidal 
suspensions. Neglect of surface effects in the computation of the colloidal virial term leads to 
spurious and paradoxical results. This pitfall is one of the several facets of the danger of a naive 
implementation of the so called One Component Model, where the micro-ionic degrees of freedom 
are integrated out to only keep in the description the mesoscopic (colloidal) degrees of freedom. 
On the other hand, due incorporation of wall induced forces dissolves the paradox brought forth in 
the naive approach, provides a consistent description, and confirms that for salt-free systems, the 
colloidal contribution to the pressure is dominated by the micro-ionic one. Much emphasis is put 
on the no salt case but the situation with added electrolyte is also discussed. 



I. INTRODUCTION 

In a complex mixture where several species with widely different characteristic time and space scales coexist, it 
is common practice to resort to a coarse grained description integrating from the partition function all degrees of 
freedom that do not belong to the main (larger) constituent [J [3 S 0] ■ This leads to a state dependent effective 
Hamiltonian for the main constituent, thereby allowing a One Component Model (OCM) description. The motivation 
for such a procedure is not only to facilitate contact with experiments, where most of the time the small constituents 
cannot be probed directly, but also to simplify the theoretical treatment. Indeed, one can then use the well developed 
statistical mechanics tools from the theory of simple liquids to study the OCM. This transposition from simple to 
complex fluids is however paved with practical difhculties, see e.g. [s, It is the purpose of the present paper to 
discuss one such pitfall arising in the context of charged colloidal suspensions. 

The system we will consider is made up of Nc charged spherical hard particles (colloids) immersed in a solvent 
with dielectric constant e, which fills a box with volume V limited by a neutral hard wall. The colloid's interior is 
assumed to have the same dielectric constant as the solvent. Each colloid bears a charge ZcC where e is the elementary 
charge and 3> 1. The medium outside the container is a structureless dielectric continuum with dielectric constant 
e' . To ensure electroneutrality, the solution contains NcZc microscopic counterions, assumed monovalent. Additional 
microions may also be present due to the dissociation of an added salt and the total number of microions is denoted 
-^micro- The particlcs interact through Coulomb forces and hard sphere exclusion, although in the subsequent analysis, 
the hard-core interaction will turn out to be irrelevant. 

The paper is organized as follows. In the situation where e = e', we consider in section |lT] the different virial 
contributions to the equation of state. In the salt-free case, we argue that the colloidal contribution Poem has to be 
negligible compared to the microionic one (Pmicro)- We then show that a naive implementation of the OCM picture 
leads to a violation of the constraint Poem <C Pmiero- Sections IIIII for e = e' and IIVI for e' < e are devoted to the 
resolution of this apparent paradox. It will be shown that in a closed cell, the surface contribution to the colloidal 
virial Poem is comparable to the bulk term, while only the latter is considered in the naive picture. Hence its failure, 
resulting from a gross overestimation of Poem- As a consequence, the knowledge of a good effective potential in the 
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bulk is insufficient when it comes to directly computing the colloidal virial in a closed box. Concluding remarks are 
drawn in section |Vl where we discuss in particular how the effective potential can be used -indirectly but from a 
standard procedure- to compute the total pressure of the system. While most of the analysis holds without salt, the 
situation of an added electrolyte is also briefly addressed. 



II. EQUATION OF STATE, EFFECTIVE INTERACTIONS AND ONE COMPONENT MODEL VIRIAL 



A. The equation of state 



We start by the simplest situation where e — e' and consider all charged species in the solution. The virial theorem 
allows to write the total osmotic pressure P (with respect to pure solvent) in the form 

where — \/{kT) is the inverse temperature and the summation runs over colloids and microscopic ions, therefore 
involving Nc + A^micro terms. In Eq. ([1]) the angular brackets denote a statistical average (that coincides with 
time average) and Fj™* is the (internal) force exerted on particle i at position r^, due to hard core and Coulombic 
interactions in the solution. Adding the force exerted by the wall to F™^ would therefore provide the total force Fl°^ 
felt by particle i. We note here that it is possible to express the pressure in Eq. ([1]) as a surface integral over the wall 
of the total (colloid + micro-ions) concentration. Applying the virial theorem to the microions only, we have 

^microfcT+j/ Y. n--F;*°M = = 7V,„ierofcT+i/ ^ r^.p'^A-^/j p^,„,{r) r ■ s\ , (2) 

XL^inicro I \i^raicro I ^ hovL I 

where the surface integral with normal oriented outward runs over the box confining the system. Inserting the latter 
equality into ([1]), we obtain 

P = p,kT + Poem + Pmicro with Poem = ( J" ^ ' i^'"* ) ] -Pmiero = ^ ( <f /Omiero(r) r ■ s) , (3) 
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where pc = Nc/V and Pmiero(^) denotes the total microion density at point r. Within mean-field approximation, this 
equation may be found in il| . The first term on the right hand side of Eq. ([3]) is the colloid ideal gas term which 
can be safely neglected in practice for the parameter range of interest here (see below) . The second term -of central 
interest here- is the colloid-colloid virial contribution and is indexed by the subscript OCM since it would be the only 
term considered (apart from the ideal gas one) in the OCM approach, restricted to the mesoscopic degrees of freedom 
{fi}i<i<Nc- Indeed, the statistical average (. . .) may be performed in two steps : 





col-{-niicro \ ' col 




(4) 



where we have introduced the microion averaged effective force F^^ exerted on colloid i for a given colloid configuration. 

The third term in ([3]), Pmiero, accounts for the direct coupling between colloids and microions. In principle, this 
third term is to be averaged over the colloidal degrees of freedom. However, even at simplified or mean-field level, a 
full iVc-coUoid simulation is computationally demanding 0, [1,0) and further simplifications are helpful. Of particular 
interest are two such simplifications, both belonging to the Poisson-Boltzmann family, that reduce the initial A^c-body 
problem onto a Nc — 1 body situation. The first one is the common cell model approach originating from a solid 
state point of view where the Wigner-Seitz cell around a colloid is constructed and then "sphericalized" for the sake of 
simplicity. The Poisson-Boltzmann equation is solved within this cell, and from the microionic density profile one can 
then estimate Pmiero- The second model is the renormalized jellium model [lo| where a liquid state point of view is 
adopted: the colloid-colloid pair distribution function (?cc(f ) is considered structureless so that other colloids around a 
tagged macroion behave as a continuous background. The charge of this background is a priori unknown, and enforced 
to coincide with the effective charge. This self-consistency requirement leads to a unique and well defined effective 
charge It has been shown that for salt-free suspensions, these two models - cell and jellium - both lead to a 

pressure Pmiero that is in excellent agreement with existing experimental data [Tlj and primitive model simulations 
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for P [Hjlill, see e.g. [13, [13, [iBl • We note that Pmicro may be coined a "volume" term 0, since - at least, 

within the cell model and jellium approaches - it does not depend on the colloidal degrees of freedom but only on the 
mean colloidal density. The good agreement one obtains with the exact pressure P for both models implies that for 
salt- free systems P ~ Pmicro- This is corroborated by a recent study of finite stiff-chain polyelectrolytes fl9'|. From 
Eq. ([3]) where the ideal gas contribution {p^kT) is neglected, this may be transposed into the following requirement: 

Poem -^micro- (5) 

A similar conclusion was reached in Ref . [20i ■ 



B. Effective interactions 



Both Poisson-Boltzmann cell and jellium approaches are not only useful to estimate the pressure, but also to derive 
effective parameters for solvent -|- microions averaged colloid/colloid interactions. By construction, the effective 
potential is that which leads to the correct colloid-colloid pair structure encoded in the potential of mean force Qcc, 
assuming pair-wise colloid- colloid interactions within the OCM model (see e.g. Although the effective potential 
has a clear-cut definition, there is no rigorous operational route to construct this object. In general, when microionic 
correlations do not invalidate the mean-field picture p^ . a good approximation is to write the effective potential as 
a sum of pair-wise Yukawa terms of the form 

o t\ 72 X /exp(K,,,.a)y exp(-K,„.r) 

with a the colloid radius, Xb = /Je^/e the Bjerrum l eng th, and Z^^ and the effective charge and inverse screening 
length computed within the cell or jellium model [U, HI, IHI- Such a "DLVO"-like expression [l|, 0, 0] would 
accurately reproduce the large distance interaction of two colloids in a salt sea P, H, 01 • Its relevance in the no-salt 
case will not be discussed. As will become clear below, we are interested here in orders of magnitude, that should not 
depend on the precise form of ([6]). 



C. An apparent paradox 

Within the jellium model, the salt-free equation of state takes a particularly simple form 

/3/micro — Z^n Pc- (7) 

Within the cell model, this expression is not exact but approximately correct. For a highly charged macroion, one has 
Z^tf ^ 1 which allows to neglect the ideal gas term in ([3]). In spite of its simplicity, the expression /SPmicro = Pc Z^.tt 
hides a complex density dependence through Z^a and is in excellent agreement with the exact pressure P found 
experimentally or in primitive model simulations, as emphasized above. In addition, the effective screening length 
reads [13] 

K^ff = AttXe Pc ^ctt (8) 

The constraint embodied in Eq. ^ may therefore be rewritten 

/3Pocm < ^cffPc (9) 

Alternatively, in the low electrostatic coupling regime (where Z^^ coincides with Zc), one should recover the ideal gas 
pressure PP ~ pc{l -\- Z^^t). Given that in this limit, /JPmicro — PcZ^a, we recover the requirement that will be 
an important benchmark for the following analysis. We now turn to the formulation of the apparent paradox. 
In the bulk of the suspension, the effective potential ([6]) provides the effective force acting on a colloid i 



(10) 



r—Vi — r 7 



Considering naively that Poem appearing in ^ and ^ is dominated in a very large system by its bulk behaviour, we 
insert pU)) into @ to approximate Poem by P*^^ with 

K.m = ^.ly.n-F^r) =jyif:^^yF^i) ' (11) 

7 = 1 / col 
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where = Ti — rj . We will subsequently omit the subscript "col" indicating the degrees of freedom involved in the 
average. Introducing the colloid-colloid pair correlation function (/cc('')i we can write 

PP:...--4 r 9Ur)'-^rd^r (12) 



6 Jr=2a dr 

= ^"^gfi^^ |l+ j + l r [5..(r)-l](l + ..„..)/^Wr)rfV (13) 

l^ia I 3(l + Krffa)^J 6 Jr=2a 

To estimate the above quantity, it is sufficient to keep the dominant term only, which is the first one on the rhs, 
arising from the long-range behavior of the pair correlation function {g^c ^ 1 at large distances). In this term, the 
curly brackets may be safely approximated by 1 since at low densities, Keft^ 1- Remembering Eq. ([5]), we obtain 

p* ~ \-L £« £ 

ocm — 2 V / 

- ^^offPc, (15) 

The factor 1/2 which appears is classical (see e.g. [H). The important point here is that estimation (fT5|) by far violates 
the constraint ([HI). A similar conclusion would be reached including the first correction in Z^^^ e:}cp{—K^iir)/r to the long 
distance behaviour g = 1 when computing the integral on the rhs of (fT^ : this yields -Poem — -^otfPc/2[l -I- ©(Keff A_b)] 
with KetfAs <C 1 in the dilute limit. The paradox here is that the very same approach that provides a contribution 
Pmicro Very close to the total pressure, gives an effective potential that apparently spoils the previous agreement, by 
grossly overestimating the colloidal virial contribution to the pressure. We will see that this feature is not ascribable 
to a failure of the functional form of Eq. ©, which provides a decent approximation for the quantity P^cm- 

D. How can the paradox be resolved ? 

The root of the paradox reported above is that approximating Poem by Poem incorrect, while Po*cni provides 
a reasonable estimate for the bulk contribution to Poem, surface effects make that in the vicinity of the wall, the 
effective force felt by a colloid differs from (fTU]) . These surface induced terms play a key role here and contribute 
a large amount to the colloidal virial Poem, no matter how large the system is. It turns out that bulk and surface 
induced contributions almost cancel each other, so that the resulting expression for Poem is much smaller than P*^^ 
and therefore fulfills the requirement ([9]). Our goal in the remainder is to illustrate this cancellation explicitly, from 
a correct description of confinement effects. To this aim, it is judicious to simplify the problem by considering the 
limit of point colloids (a — 0), and by identifying the effective charge with the bare one Zc- Considering charge 
renormalization effects is here immaterial and focussing on dilute systems where ko is small, finite a effects do not 
affect our main conclusions. In the bulk of the suspension, the effective potential therefore takes a simple Yukawa 
form 

f3v,sir) = ^2^^ exp(-^r) ^ ^^^^ 
r 

with = AttXbZcPc- 

At this point, a comparison with simple electrolytes seems appropriate, for the aforementioned cancellation is 
already present. For our discussion, we may consider that the role of the colloids is played by the cations, and that 
the anions constitute the remaining "microions" . The pressure has to be close to [ll, |j] 

^3 

/^Pclcctrolytc — Panion ^ Pcation TT^j 7 

(17) 

z47r 

with equal mean densities Panion — Pcation- From the contact theorem, we deduce the densities at the wall 

^3 

Pa.nion{wall) = Pca.tion{wall) ~ Panion - T^— . (18) 

4o7r 

Rewriting ([3]) in the form 

/5Polcctrolyte = Peations + ^ ( ^ n • i^- "* ) + Panion (wa^O (19) 
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we obtain from (flTl) and ITS' 



E 

\i^cation 



int 



K 



(20) 



Given that = SirXsZ^pcation, we have K'^/{PP*^^^) oc kXb which is a small quantity for a dilute system. We 
explicitly see here that the "colloidal" virial [Ihs of (|20p up to a factor (3] is by far smaller than the estimation P*^^- 



III. WALL MEDIATED FORCES WITHOUT DIELECTRIC DISCONTINUITY 

In the vicinity of the wall, the colloids do not see a spherically symmetric environment. As a consequence, 

1. the usual exp(— Kr)/r pair interaction is modified. 

2. the mean force acting on a colloid does not vanish. This is a one body, wall induced effect, mediated by the 
microions. It is therefore an internal force, that should be taken into account in It should not be confused 
with the external (and short range) direct colloid-wall interaction. 

Evaluating the rhs of ([4]) therefore requires a careful computation of both types of microion averaged colloidal forces. 
To this end, we need the solution (j>z{p, z') of Debye-Hiickel equation V^02 = k'^0{z')4>z in the case where a test charge 
is located in the solution a distance z from an infinite neutral wall. We have introduced the Heaviside function 9 
and cylindrical coordinates (p, z') such that the test particle is located at (0, z) with z > 0. The planar geometry 
approximation for the wall is sufficient provided the cell size or radius of curvature is much larger than Debye length 
1/k. We start by the situation of equal dielectric constants inside and outside the solution (e = e'). The electrostatic 
potential rnay be written in the form of a Hankel (two dimensional Fourier) transform [26)] where q and p are conjugate 
quantities f24. [25| 

0,(p,z') = Z,XbJ^°^ (Irf ''''''''''^ ^ e~''\'-''\^^Joiqp)qdq ; k = ^^^T^ (21) 

The second term in the integrand (e^''!^^^ I) gives exactly ZcXb Gxp{—Kr)/r where r — [p'^ + {z — z')^]^/^ is the 
distance to the source. This is the standard Debye-Hiickel potential which dominates in the bulk. The remaining 
term, which vanishes at large distances {kz or kz' ^ 1) is due to the presence of the interface. 



A. One colloid ion average force 

The force felt by a colloid located a distance z from the planar interface follows from ([2T|) , considering the electro- 
static potential (jfz = (f'z — Z^^b e~'^^ /r where the self term has been subtracted: 



^F^.^aii = nZ,^/?04O,z') 



= ^c^^ I ^e-^"^ qdq i-n). (22) 
k + q 



In this equation, n denotes the unit vector perpendicular to the interface pointing outside the solution. We coin the 
force "colloid-wall" and for notational convenience, we henceforth omit the superscripts "int" and "eff" . This 
force repels the colloid from the wall [k = (k^ + g^)^/^ > q], as a result of microions imbalance between the half of 
the colloid exposed to the wall, and the other hemisphere. Inserting (|22p into ^ we have 

^ (E^^'-^*-™'^") = ^ / / Pciz)r-F,^^,aii{z)dz. (23) 



■wall 



To leading order, the above integral may be computed assuming a uniform density of colloids Pc{z) — Pc- In (j23p . 
r denotes the absolute position with r = s zn (s is therefore the orthogonal projection of r onto the wall). We 
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neglect the term in — zii (that would contribute proportionally to the surface of the system), so that 

37 (^Jrm-n,ai}l - ^ f j s-F,^^aii{z)dz (24) 



II Jo 



^^^^I^bH s-nd^s) Td. r^l^^e-^^V^.d, (25) 
\Jwall ) h h ^/W+q^+q 

~ -\p,ZIkXb (26) 

D 

Incidentally, this is exactly the Debye-Hiickel form for the excess pressure of an electrolyte [see Eq. P7|) ]. For dilute 
systems, this quantity is small compared to PcZc^ as emphasized earlier. The constraint Q is therefore fulfilled. 

B. Colloid-colloid interactions 

Within the simple Debye-Hiickel treatment, the potential of interaction between two colloids near the wall (one at 
2, the other at z', with a lateral distance p between them) is Zc4>z{p, z') = Zc4>z'{p, z). To calculate the force felt by 
the colloid at z due to all neighbors, we assume again a uniform distribution of neighbors: 



f3Fcoi-coi{z) = npc dz 2tt p dp - 

Jz'=0 JO 



dz 



(28) 



= ~np,Z',XB / dq / 27:pdp ^^_-L ^ - 1 e'^-V^^+r —== Mqp)qdq (29) 
Jo Jo \ V + q^ + q J V k2 + q2 

The component of the force parallel to the wall vanishes upon averaging. 

Inserting this force into ^ and proceeding along similar lines as in Eqs. (|24p sq, we have 

4j (yrrF,-cou) ^ P^Z^^Xb H 2npdp T (- ^£±^'^ + 1^ Jo(gp)^dq. (30) 



\^ / Jo Jo \ V'* +q^ + q J 1^ +r 

Both expressions (1^^ and ([50)1 are of the form of a Hankel transform at the origin q = of the inverse Hankel 
transform of a function A{q), with A = {. . . - l)e^''^^/k in (HH) and A ^ {- . . . + in ^0\j. This is nothing but 

A(0) 26], which vanishes in both cases. Therefore, with the approximations proposed, the force in (|29p and the virial 
term in (j30p vanish. To be more specific, we compute explicitly the integrals in pop : 

^ (Y.^^■F^-col}| ^ ipc^c(-l + l). (31) 

The term in +1 in the parenthesis arises from the term in +1 in Eq. ([30| . which gives the usual "bulk" e^"'/?' 
pair interaction, as already mentioned. The associated virial is /3-Po*cm = as obtained in (fTSll . The present 

calculation shows that this term is canceled by an opposite wall induced contribution. If the simplifying assumption 
gcc = 1 is relaxed, the resulting expression for (|31[) no longer vanishes but remains negligible with respect to PcZc- 
On the other hand, relaxing the assumption of a uniform profile Pc{z) leaves the result unaffected, as will be seen in 
section |lVl 

We conclude here that summing the two contributions from Eqs. (^7)) and (Plj) provides a value for Poem that is 
compatible with the constraint ([9]). 



IV. ANALYSIS IN PRESENCE OF A DIELECTRIC DISCONTINUITY 



In this section, we extend the previous analysis to the situation where the dielectric constants are not matched: 
T] = e'/e ^ 1. The relevant parameter range corresponds to 77 < 1 e.g for water droplets in air in a spray-drying 
experiments. The first important difference with the ry = 1 case is that the equation of state ([3]) takes a different 
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form. The pressure is indeed not solely given by the contact densities of charged species at the wall, but contains 
additional electric contributions (polarization or image effects). On the other hand, Eqs ([T|) and ^ are still formally 
correct provided one also includes in the "internal" forces the electric forces from the wall. The resulting equation of 
state reads 



Here 



T"' = — E^I - —E®E (33) 
Stt 47r ^ ' 



is the Maxwell tensor, with E the local electric field and I the isotropic tensor. 
The counterpart of ([2T|) now reads: 



<PM^') = Z,\b r ('^^e-'=(-+^') + e-^\^-^'\]\Mqp)qdq ; fc = T^^T^. (34) 
Jo \k + m J k 

As in the case rj — 1 and as long as ?7 7^ 0, the corresponding interaction between two colloids decays as p^^ at large 
distances parallel to the wall (see [13] for a discussion of this dipolar- like term). When 77 = 0, the wall can be formally 
removed considering the electric image located symmetric to the z — Q plane. 

The colloid-colloid and colloid- wall interactions readily follow from (|34p. At short distances z ^ 0, the latter 
diverges like z~^(l — Ti)/{\ -I- 77) [25[ . which corresponds to the unscreened interaction of a particle with its own 
image. This divergence means that the uniform colloid density cannot be invoked when it comes to computing (j23p . 
To obtain the leading order behaviour, we can assume that the colloids are distributed with the Boltzmann weight 
Pc{z) = pce-yiY>[- I34>c-wan{z)], where Fc-waii = -'^4>c-waU and the potential (j)c-waU deriving from ([M]) vanishes for 
z ^ 00. The precise knowledge of this potential is however not required since 

^ ( V'n'-Fi-w.ai; ) = [ d^S I r-Fc-n,aii{z) eyi^[-f3(t)c-waii{z)]dz. (35) 

I '^^ J wall Jo 

~ PckT [exp(-/30c-tt>aH(2))](^ (36) 
-PckT. (37) 



This term therefore cancels the ideal gas one on the rhs of ([321) . 

The wall induced colloid-colloid contribution to the colloidal virial may be computed along similar lines as in section 
IIIIBI An expression involving again a Hankel transform composed with its inverse is again obtained, with now a 
function 



A{q) = ^^I^-^^ 1 Hdz pc{z)e-^^ (38) 

X^ti^ + q^+riq I \J + q^ Jo 



Since A{0) = 0, we conclude here that 



) ~ 0, (39) 



so that the total colloidal virial [including colloid-colloid and colloid-wall interactions] is close to —pckT, which is a 
small quantity compared to the microionic contribution ZcPckT. Equation p2p can finally be rewritten 



V. CONCLUDING REMARKS 



Before briefly discussing the situation where a salt is added, two comments are in order. 
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A. Closed cells versus periodic boundary conditions 

From the previous discussion, it appears that the equation of state ([3|) holds when the system is confined by a hard 
wah, and would fail if periodic boundary conditions (pbc) would be enforced. The inadequacy of P*^^ to approximate 
Poem may then be phrased in the following way 



3VP* 

^ * ocm — 2 



E^^^-nV..-i2n)) (41) 



= l " ' pbc 

+ {^v^-Ff\ (42) 

\ ^ / hard walls 

where in (j4ip . the sum involves all periodic images of the cell considered: n is a vector with components in Z,^, which 
indexes the center R„ of a given image of the "central" cell. The central cell has Rq = and since we deal here with 
a short range effective potential, the sum over n may be truncated to retain only the 7 terms with |n| < 1. 
However, for any simple fluid where the forces Fi are given, (j42p would be an equality. Indeed we have 

(\ simple fluid i v simple fluid 

T.^rF^) (43) 
« / hard walls \ *>i / hard walls 

where the rhs shows negligible dependence on the boundary conditions provided the system is large enough, and can 
then be computed with pbc provided the correct forces are considered [Fi — Fij{rij — Rn)]- Hence 

(V simple fluid i v simple fluid 

E ^» ■ ^0 =UJ2J2''^^- (^^^ - ) (44) 
i I hard walls \ *,i " / pbc 

The difference between equations (H^ and (|44p illustrates the important role of microions. We may also consider that 
the ^ sign in (|42p arises from the density dependence of the effective pair potential. 

A natural question at this point is : does the knowledge of the "bulk" effective potential (O between colloids allow 
to compute their virial Poem as it appears in ^ ? The answer is positive in a closed cell, at the OCM level, provided 
that due account is taken for the dielectric images of the colloids. In the following section, we address a related 
question, and discuss how the full pressure of the colloidal system may be recovered, assuming again that the only 
information at hand is that of the bulk effective colloid-colloid interaction. 



B. Back to the DLVO potential 

We consider here a simple liquid that interacts with a pair- wise potential given by Eq. ([S]), with effective parameters 
3> 1 and K*ff = A'k\bZ*^^p*^ (salt-free case, for simplicity). These parameters are fixed a priori, and chosen to 
coincide with those relevant for a colloidal suspension at pc — p^- The potential of interaction is therefore density 
independent and the system, later referred to as "auxiliary" , can be studied for pc =/= p*. 

We consider the parameter range (essentially low density) where the excess pressure of such a system is well 
approximated by Po*cm i^i Eq. ([14]): 

^p* _ 27rpgZ*f Aj3 _ 1 pI 

'^ate ^ Pc 

Incidentally, the contact theorem indicates that the contact density in the case where the system is confined by a 
closed box, reads pciwall) ~ pl^Z*^^ / {2p*^). This quantity is much larger than the mean density pc (except when pc is 
extremely small, a limit of little interest here) . This excess with respect to the mean density is to be contrasted with 
the depletion from the wall that is present in the original colloidal system containing microions: Eq. (|22p for e — e' 
shows a repulsive colloid-wall behaviour, and the depletion is even stronger when e' < e due to like-sign images, see 
the discussion after Eq. (p4|) . 

The pressure of the simple liquid with DLVO interactions, close to Po*cm: l^as a priori nothing to do with the pressure 
porigmai ^j^g ^.gg^j coUoidal systcm. It has also nothing to do with the colloid virial contribution entering Eq. ([3]). 
However, for Pc = P*ci the colloid-colloid structural information is the same for both original and auxiliary systems. 



9 



One may then invoke Kirkwood-BufT identity '28'| which states that the inverse compressibiUty of the original colloidal 
suspension coincides with the long wave-length limit of the colloid-colloid structure factor Scc{k) '■ 



X 



dpc 



ScciO). 



(46) 



The compressibility in our auxiliary simple liquid with fixed potential of interaction is therefore the same at pc = p* 
(and only at this density) 



^^original 



dpc 



dPZ 



dpc 



(47) 



This offers a means to compute the equation of state of the original colloidal system from integrating the inverse 
compressibility of the auxiliary one. In this integration, due account must be taken of the density dependence of 
both and . The previous integration procedure therefore requires to consider the auxiliary system for several 
values of Pc for a given p* [to compute the derivative in the rhs of (|47p ] , before scanning the range of interest for p* . 
Of course, the general procedure outlined here does not depend on the specific form of the effective potential, and is 
equally valid when salt is added. It turns however that the DLVO potential together with the salt-free approximation 
([14]) -which leads to (|45)) - provide a clear illustration of the procedure. From (|45l) . we obtain the rhs of Eq. (|47| : 



d(3P* 



dpc 



at Pc = Pc 



(48) 



To compute the Ihs of (|Tf|) . we may come back to the jellium model which gives /3poi'g'"ai ^ Z^^pc- In this expression, 
the effective charge may depend on the density, but for salt-free cases, this dependence is at most logarithmic for 
Pc — *■ [lol | and provides only a subdominant term to the compressibility, so that 



dl3P 



•original 



dpc 



Evaluating this expression at pc — p* where Z^a — Z*^^ , we recover Eq. 
(jiT]) but also the consistency of the underlying DLVO potential. 



(49) 

48|) . This not only illustrates the identity 



C. Situation with added salt 



When the suspension is dialyzed against a salt reservoir, most of the technical analysis carried out earlier is still valid. 
We consider a similar auxiliary system as in section IVBl with effective screening length such that k*J > 47rAB2'*jj.pc 
due to the screening by salt ions ^2^. The effective charge and screening lengths are again chosen to coincide with 
those of a colloidal system at a particular density p*, but are otherwise density independent. Equation (|47p still holds 
while Poem is given by Neglecting again the integral on the rhs of ([T^ . and inserting the resulting Poem (33 1 
we obtain: 

original 



where we have replaced Z*^^ by Z^ff and k,*^^ by ^^^ff after computing the rhs of ([47)l . Here, the prefactor A reads 

■4 = 1+ J":""^' (51) 

3(1 -I- Keff a) 

Is relation §0^ compatible with Poem < P = p°"ginai 7 Neglecting Poem (together with pckT) in we have 
P ~ Pmicro which in the jellium model is given by k^^^ / (AttXb)- With the help of |23], we arrive at 



5Pm 



dpc 



^nXBPcZ]^_ (52) 



T off- 
Equations ([50]) and ([52]) give the same result provided A is close to unity, which means Kefffl < 1. We conclude here 
that omitting the colloidal contribution to the pressure. Poem, is inconsistent when K^aa > 1. It turns out however 
that A increases very mildly with KeffO (e.g it is close to 1.2 for Keff a — 4). A more precise discussion would require 
to consider the full rhs in (llSp . which is beyond the scope of this paper. Finally, we note that in the salt-free case 
where n^na = irjcZ^ifXg / A with 77^ = Airpca^ /2t the colloidal volume fraction and Z^ffX-Qja on the order of 10 for 
highly charged colloids, we have Kefffl < 3 and therefore A close to 1 even for packing fractions as high as 10%. 
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D. Summary 

We have seen that for a salt-free cohoidal suspension, the coUoidal contribution Poem to the equation of state [as 
written in Eq. ([3])] is a neghgible quantity. This feature may easily be overlooked in a naive implementation of 
the One Component Model, where only -Poem; the bulk contribution to Pocmj is computed. The fact that P*^^ is 
of the same order of magnitude as the total pressure P of the suspension, is not compatible with the requirement 
^ocm ^ ^micro — that has emerged as a central constraint in our analysis. We have shown that no matter how 
large the system is, surface effects that require the resolution of Poisson's equation in the vicinity of a confining wall, 
contribute a large amount to Poem- To zeroth approximation, these surface terms cancel the bulk value P^cm^ so that 
one finally recovers Poem <C P. 
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